一种超导重力仪磁悬浮力的等效计算

为实现超导重力仪磁悬浮力的精确计算,以GWR型超导重力仪为模型基础,采用有限元的思想,将超导球表面电流理想化为多个等高共轴电流环,计算出各个电流环与超导线圈的作用力,求和得到线圈与超导球间的磁悬浮力。利用MATLAB完成计算程序实现,通过改变下线圈电流和上、下线圈电流比,获得满足一定条件的磁悬浮力及其梯度。选取合适的模型参数,计算出线圈对质量为m=4.069 g超导球的磁悬浮力大小为:Ftotal=3.988×10^-2N,磁悬浮力梯度为:-9.699×10^-3N/m,此时悬浮力梯度合适,满足系统稳定性和灵敏度的要求。     

中国聚变工程实验堆真空室超压保护系统中爆破片的选型与计算

参照相关标准，对聚变装置真空室超压保护系统(VVPSS)中爆破片进行了选型。结合VVPSS的工作要求，完成了爆破片的设计计算，初步得到爆破片直径为882mm，厚度为1mm。利用有限元分析软件对多种型号爆破片进行结构分析比较，最终选用了反拱环向开缝型爆破片。对最终选定的爆破片进行优化设计，使其达到设计要求。     

高温、高压、高应变速率动态过程晶体塑性有限元理论模型及其应用

对于高温、高压、高应变速率加载条件下的材料冲击变形行为,动态晶体塑性模型能够直接反映晶体中塑性滑移的各向异性及其对温度、压力和微观组织结构的依赖性,因而广泛应用于材料的动态冲击力学响应、微观结构演化以及动态损伤破坏的模拟。本文综述了高压冲击下动态晶体塑性有限元的理论模型,主要包括变形运动学、包含状态方程的超弹性本构模型和晶体塑性本构模型,涉及位错滑移、相变、孪生等塑性变形机制,以及层裂、绝热剪切带等动态破坏方式。     

非定常Navier-Stokes方程基于两重网格离散的有限元并行算法

基于两重网格离散和区域分解技巧，提出三种求解非定常Navier-Stokes方程的有限元并行算法.算法的基本思想是在每一时间迭代步，在粗网格上采用Oseen迭代法求解非线性问题，在细网格上分别并行求解Oseen、Newton、Stokes线性问题以校正粗网格解.对于空间变量采用有限元离散，时间变量采用向后Euler格式离散.数值实验验证了算法的有效性.     

基于像差特征分析的变焦光学系统装调

提出了一种基于像差特征分析的变焦系统共轴性装调方法。针对某型20倍变焦光学系统的装调，在分析变焦光学系统装调要求的基础上，通过采用光学定心加工技术，提高变焦系统各镜组内光学元件的同轴度；采用Zygo干涉仪，检测变焦系统在长短焦位置的像差特性分布；借助CodeV软件仿真计算各光学元件在系统中的公差灵敏度分布，并确定产生像差影响的敏感光学元件；在中心偏测量仪上，完成最终光学系统像质的调整。此外，还设计了一种变焦光学系统各动组间同轴度调试装置，对变焦相机主镜筒机械内孔轴线与直线导轨的平行性进行了精确测量，保证了动组组元光轴的同轴精度。装调结果表明:变焦系统成像质量有明显改善，像面一致性得到保证，长短焦轴上传递函数值分别优于技术要求值0.45和0.55，长短焦轴外0.7视场传递函数值优于0.25和0.35，实现了高精度装调，验证了该方法的可靠性。     

Investigation on the effect of density ratio on the convergence behavior of partitioned method for fluid–structure interaction simulation

In order to investigate the effect of density ratio of fluid and solid on the convergence behavior of partitioned FSI algorithm, three strong-coupling partitioned algorithms (fixed-point method with a constant under-relaxation parameter, Aitken’s method and Quasi-Newton inverse least squares (QN-ILS) method) have been considered in the context of finite element method. We have employed the incompressible Navier–Stokes equations for a Newtonian fluid domain and the total Lagrangian formulation for a non-linear motion of solid domain. Linear-elastic (hyper-elastic) model has been employed for solid material with small (large) deformation. A pulsatile inlet-flow interacting with a 2D circular channel of linear-elastic material and a pressure wave propagation in a 3D flexible vessel have been simulated. Both linear-elastic and hyper-elastic (Mooney–Rivlin) models have been adopted for the 3D flexible vessel. From the present numerical experiments, we have found that QN-ILS outperforms the others leading to a robust convergence regardless of the density ratio for both linear-elastic and hyper-elastic models. On the other hand, the performances of the fixed-point method with a constant under-relaxation parameter and the Aitken’s method depend strongly on the density ratio, relaxation parameter selected for coupling iteration, and degree of deformation. Although the QN-ILS of this work is still slower than a monolithic method for serial computation, it has an advantage of easier parallelization due to the modularity of the partitioned FSI algorithm.     

基于改进格林元的压裂水平井动态分析方法

以建立高效的动态分析方法为出发点，以边单元作为求解点，改进传统的格林元方法，减少未知数和求解矩阵维度；并提出基于改进格林元的加密网格加密方法，保证考虑复杂裂缝网络的压裂水平井动态模拟的早期精度.退化模型与半解析解、数值模拟结果进行对比，验证本文基于加密网格的改进格林元方法的准确性和动态分析的高效性.最后进行动态响应的敏感性分析，结果表明：①格林元方法是一种高精度的动态模拟方法，将求解节点设置在网格的边上可以提高压裂水平井动态模拟的速度；②改进格林元方法的加密基于叠加原理，不需要通过插值近似，其求解精度高.在相同加密网格条件下，基于本文改进格林元方法的加密效果比有限差分加密效果更佳；③复杂裂缝导流能力、改造区渗透率提高倍数、改造区大小等参数对压裂水平井动态特征影响较大，在动态分析和参数反演时，应着重考虑这些因素的影响.